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The properties of aqueous fluid of potassium fo nn atewere studied theoretically and 
experimentally in order to find arT improved secondary refrigerant for indirect 
refrigeration systems. The most important advantages compared to t he traditional 
aqueous solutions of alcohols and glycoles are good thermodynamic properties, low 
toxicity and non-flammability. A remarkable benefit is a reduced change of a laminar 
flow. Volumetric heat is lower, which brings about a higher rale of mass flow, if the 
temperature change of the fluid is maintained constant. However, because of the clearly 
lower values of viscosity, the pumping power demand is at same level. Good thermal 
conductivity promotes good heat transfer, rr? 1997 Elsevier Science Ltd and II R. 
(Keywords: heat carrier; chemical composition; physical property; viscosity; specific heat; toxicity; 
flammabiHry) 



Le formiate de potassium comme frigorigene secondaire 

Pans le but dc trouver tut tnedleur frigurigene secondaire pour ies sys femes de froid indirect, 
on a fait I 'etude theorique et e.\ pet imen tale des proprietes dune solution aqueuse de 
formiate de potassium. De bonnes prapr ivies thermodynanuqites. unc faibfe to.xitite et 
rinflammabiliie sont ses principally attntts par rapport att.x solutions aqueuses tradition- 
nclles d'alcools vt de g!ytof\ Ln faeilite d'eviter I'ecnulement taminaire est un a vantage • 
remarquable . So moindrc eludeur volumetrique necessite un debit massique super ieut pour 
unc difference de temperature du ftuide donnee. C'ependant, en raison de so viscosite 
nettentent infer ieure . la putvstnn e au pompagc reste au inctne niveau. So bonne conduclivitc 
thertnique favorise un btm tnm\fert de chalcur. ({': Elsevier Science Ltd et II R. 

(Mots des: caloporlcur. composition <Jiimit|tte; potassium: proprietc physique, viscositc; chalcur 
massique: loxicile; iiiflamtn.ihiiiu) 



Introduction 

The requirement to minimize (he use and release of 
hamifuj^-gr^ hazardous refrigerants leads more and 
< m"6^fe^o"com TpTi cTTel rTfCttrtterTTp a cfca ye sTT rope r t i e s of 
heat transfer fluids or secondary refrigerants, as t hey 
are often called, arc essential for a reliable indirect 
refrigeration plant. Besides thermal and transport 
properties, other important properties arc cnrrosivii v. 
safely. toxicity. environmental characteristics 
and compatibility with common construction 
materials. This paper introduces a mm traditional 
suhstancc. potassium formate (IT), ami -.-iMiiparc* 



its properties with some most frequently used aqueous 
fluids. 

A number of papers on alternative secondary fluids 
have been published in the last few years (see e.g. Refs 
I and 2). Also some new fluids have been proposed, 
such as rape oil 2 as an environmentally friendly 
alternative for soil heat pump systems. The majority 
of traditional fluids arc based on aqueous solutions of 
salts or some organic liquids, alcohols and glycols 
being the most common choices. Some non-aqucous 
liquid*; are used in low temperature applications. 

None of the alternatives available can fulfil all the 
requirement* for a good fluid Aqueous salt solutions. 
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Nomenclature 






A 


Heal transfer area 


<\ 


Heat transfer coefficient 


c 


Specific heat 


A 


Thermal conductivity 


tt 


Diameter 


<t> 


Heat flow 


G 


Thermal conductance 




Efficiency of the pump 


L 


Length 


V 


Kinematic viscosity 


th 


Mass Maw rate 


P 


Density 


Nu 


Nussclt number 


^ 


Friction factor o 


P 


Pumping power 






Ap 


Pressure loss 






P 


Pressure 


Subscript. 


r 


Pr 


Prandll number 


0 


Reference state 


R 


Thermal resistance 


a 


Secondary refrigerant in (he air cooling 


Re 


Reynolds number 




coil 


U 


Overall heat transfer coefficient 


e 


Secondary refrigerant in the evaporator 


w 


Flow velocity 


r 


Secondary refrigerant 



like sodium chloride or potassium c hloride, are very 
corrosive and a carcful^seleciion oT inhibitors and 
materials is necessary. In practical installations, 
however, corrosion may proceed faster than expected 
and a frequent monitoring of the system is necessary. 

The solutions of the most usual organic compounds 
have low thermal conductivity and high viscosity, 
especially at low temperatures. Therefore, their heat 
transfer is weak and the pressure loss is high, 
particularly at lower temperatures. One remarkable 
disadvantage is a risk of the laminar flow in the heat 
exchangers and therefore a dramatic drop in heat 
transfer. This poses a particular problem to shell and 
tube heat exchangers and plate finned air cooling heat 
exchangers. To increase the mixing of the flow, 
twisted tape inserts are often used in the lubes, but 
this causes further pressure loss. A plate type heal 
exchanger is not so sensitive because of a turbulence 
promoting form of the flow channels. Another reason 
is the low value of the critical Reynolds number, 100- 
400 compared to 2000-2500 in pipes. This transition 
from a turbulent flow to a laminar flow may be 
particularly problematic in cases where temperature is 
changing during the operation. 

In the following, the main features of the aqueous 
solution of PF are described and some results of the 
comparisons with oiher organic aqueous solutions 
given. The aqueous mixtures to be compared are ethyl 
alcohol (EA), ethylene glycol (EG) and propylene 
glycol (PG). 



Chemical nature 



Potassium formate HCOOK_is a salt olUbrmic)acid.J t 
can be prepared byabsorbing carbon mtrftbxide'in 
50-80% KOH at I00-200°C at a partial pressure of 
CO> 690kPa. The reaction is fairly slow 3 . The 
colourless and highly deliquescent crystals of PF 
have an orthorhombic symmetry and a density of 
I9l0kgm"\ The crystals melt at I67T and 



decompose almost entirely to the oxalate at about 
360 W C with the evolution of H 2 . 

Potassium formate is very soluble in water. The 
saturated aqueous solution ail8°C contains 76.8% 
PF by weight. Aqueous solution of PF may be readily 
manufactured from formic acid and KOH. 

Th e corrosive pro perties of the solution are wea^at 
the lower temperature levels. In aqueous solutions iTis^ 
weakly alkaline and thus corrosion retarding. The 
corrosion can be prevented by proper inhibitors. 
Some traditional compounds are not soluble, examples 
of the suitable inhibitors are I. 2, 4-triazole and 
borax. 

The compatible elasiomeres are. for example, nitril \ 
rubber, chlorobuiadienc rubber, chlorobutadiene )} 
rubber, polyvinylchloride, polyamidc, polyethylene f 
and aramides. 

One remarkable advantage is it is environmentally 
friendly. Potassium formate is already widely in use in 
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CONCENTRATION, w-% 

Kigure I Kjrc/ing puini tempera lure, o f fluids compa red. 
EA. EihylalcoKo"t~Ftr ethylene glycol; ^(J. propylene 
glycol; PF. potassium formate 

Figure I Temperatures Jts points de congelation de: EA. 
ethanol; EG. ethylene-glycol; PG, propylene-glycol; PF. 
jormiute de potassium 



278 



A. Aittomaki and A. Lahii 



ihe oil-drilling industry and as an ice removing agenl 
for aeroplanes. 

Thermodynamics properties 

figure I shows the freezing points of aqueous 
solutions of the fluids compared. The earlier values 
of freezing point data 4 of PF solutions were verified 
with measurements in a cryogenic bath. It was found 
that the eutectic point may be even as low as -70°C 
with a saturated solution of PF. As seen in Figure I 
PF solutions generally have lower freezing points 
than the respective glycol solutions. 

In the following, the fluids are compared with 
respect to their refrigerant properties at two different 
concentrations: the first concentration giving the 
freezing point -I0 C C and the second -25 U C. The 
respective concentrations required for these freezing 
points arc the following: 19 and 35% for EA, 23.5 and 
41 % for EG, 25 and 44% for PG, 19 and 35% for PF. 

Volumetric specific heal for the solutions arc 
presented in Figure 2. The specific heals of PF 



solutions were estimated with a thermodynamic 
estimation program 5 . The specific heat was also 
measured with a reaction calorimeter (Met tier RCl) 
at two concentrations at various temperatures above 
room temperature. The agreement between the mea- 
sured and the estimated values was good. As seen in 
Figure 2 the volumetric specific heat of PF is relatively 
low which means higher flow rates if a specified 
temperature change of the fluid has to be maintained. 

Kinematic viscosity for the solutions as a function 
of tempera lure is plotted in Figure 3. The viscosities of 
PF solutions were measured with a capillary lube 
viscometer. Potassium formate has remarkably lower 
viscosities' lhan the other fluids. 

Methods of comparison of fluids 

The comparison of secondary refrigerants cannot be 
based on one single property at a time, because 
several thermodynamic properties influence the per- 
formance. When studying the behaviour of the 
refrigeration system, a comparison can be applied: 




•10 -5 0 5 10 15 
TEMPERATURE, °C 



20 



Figure 2 Volumetric specific heat for the solutions having 
the freezing point temperatures - 10 and -25°C 

Figure 2 Chaleur massique volume trique pour des concen- 
trations donnant des temperatures de point de congelation de 
^10 et -25°C 




-10 -5 0 5 10 15 20 
TEMPERATURE, °C 
Figure 3 Kinematic viscosity of the solutions 
Figure 3 V'iscosite cinematique des solutions 



(a) to an existing system as a retrofit fluid; or 

(b) to a new system. 

In the first case, the system components, like the 
pipes, the pump and the heat exchangers remain 
unaltered, only the secondary refrigerant and the 
mass flow arc changed. This case may become 
indispensable if, for example, the rate of heat transfer 
of the old fluid is inadequately low. 

In ihe second case, ihe whole system can be selected 
and dimensioned according to the requirements of the 
fluid. The dimensioning criteria can be based on 
different principles: 

(1) An optimization of the whole system using the 
total cost as the criierium. The sizes of the 
compressor, the air cooler and the size of the 
pipes should be included in the optimization. 
This is a rather complicated way and not very 
well defined or clear because of the great number 
of the parameters. 

(2) An optimization of the system using the total 
energy input (the compressor and the pump) as 
the criteria. When the mass flow rate of the 
secondary refrigerant is increased, the mean 
temperature of the fluid is increased and there- 
fore the evaporation temperature, as well. 
Simultaneously, the energy supplied to the 
compressor is decreased and to the pump 
increased. Of course, the optimum is influenced 
by the relative amount of the pumping power. In 
this case, the dimensions of the pipes can be 
optimized separately using the total costs of the 
pipes as the criierium. According to earlier 
studies 6 , the optimum velocity is 0.5-1 ms' 1 . 
The total costs are rather insensitive to ihe 
velocity in ihe range given above. 
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Simple comparison of thermodynamic performance 

The relevant thermodynamic properties can be 
combined in different .ways in order to describe a 
particular performance. For example, the pressure 
drop or heal transfer properties can be compared 
separately, as done by Granryd and Melinder. 
However, pressure loss and heat transfer are inter- 
related in the heat exchangers, because an increased 
flow rate gives a better heal transfer with smaller 
temperature differences across the boundary layer. 

One way to compare the pressure loss and ihe heal 
transfer simultaneously is to look at the pumping 
power, which is required to give a specified tempera- 
ture diflcrence between the fluid and the heal transfer 
surface. In order to obtain u simple approximate 
formula for illustration purposes, we can use the 
Blasius correlation for the friction factor and the 
Dittus-Boelter correlation for the heat transfer (in 
pipes) in the turbulent flow region: 

£ = QRe~ 02 



Nu - Chi a Re Pr 



(I) 



From the equations above, we obtain the pumping 
power per pipe length: 

P/L^CF^d™ (2 , 

where C Pl is a constant coefficient depending only on 
Q and C^. F P is a property coefficient: 

Fp = ^pA :j "Pr- 0 ' (3) 

In the similar way, for the laminar region we obtain 

If the temperature difference between the wall and the 
fluid, as well as the heat flux are constant, the heat 
transfer coefficient .q is constant. In this way the 
property coefficient F P indicates the differences in 
pumping power demand between the fluids, as ihown 
for turbulent flow in Figure 4. In the laminar region, 
the differences are indicated directly by the dynamic 



viscosity (Equation (4)), but usually the laminar case 
should be avoided and is therefore of minor interest. 



Properties in a brine loop 

One alternative method of the comparison is some 
kind of simulation of a whole secondary refrigerant 
loop. The influence of the fluid properties can be 
found oui by applying a simple approach. A more • 
sophisticated way is not needed to show the major 
differences between the fluids. An elementary model- 
ling was carried out for the case of a simple loop 
shown in Figure. 3. The comparisons were based on a 
reference slate denoted by subscript 0. The reference 
fluid selected is propylene glycol, the reference slate is 
specified in the Appendix. The calculations are based 
on the simple correlations and formulae shown below. 

( 1 ) Heal transfer in the tubes of the evaporator and 
air cooling coil. Diltus- 
according to Equation (I) 

(2) Pumping power 



(5) 



(3) 



where the efficiency ?/ p was assumed to be 
constant (0.7) 

Heat flow in the heat exchangers 

°\n(&T l /±T 2 ) (6) 

where &T k and AY" 2 are the temperature differences 
at both ends of the heal exchanger and G is the 
thermal conductance (= AU y where U is the overall 
heat transfer coefficient). 

The heal transfer coefficients on the air side and on 
the refrigerant side were kept constant. The change of 
the heal transfer coefficient in the air cooling coil 
compared to the reference state is calculated using 
Equation (1): 

^ (0 Ck>) Grti) Co) do) 

(7) 




-15 -10 -5 0 5 10 15 20 
TEMPERATURE, °C 
Figure 4 The property factor F ? in Equation (2) 
Figure 4 Facteur do V Equation (2) 



\ 




EVAPORATOR AIR COOLER 

Figure 5 The simple cooling loop simulated 
Figure 5 Simulation de la boucte de refroidissement simple 
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For the heat transfer rale in the coil three formulae 
can be written: 



Air side: 

Liquid side: 

=m<i't(T t y - T f2 ) 
Between the flows: 

- ^r3 - ?*i2 + ^r2 



In 



7-.i - r ri 



Where the '.hermal conductance is: 
G. = AJ(R t + I /«,,) 



(H) 
(9) 

(10) 
(II) 



/?, is the (hermal resistance in the brine side from the 
wall to air. 

Similar equations can be written for the evaporator: 

(12) 
(13) 
(14) 

(15) 



* c =r" c c,(T, t - r n ) 



T„ - T. 



ln r tl -r c 

G c = AJ{R* + l/a le ) 



The change of pumping power in regard to the 
reference state can be based on Equations (I) and (5): 



V'w. \«w VP*/ 



p=p a {^\ I-) l-^l ("6) 

Only the pressure loss due to friction in the pipes is 



12000 




8000 
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Figure 6 Calculated capacity of the evaporator. The 
freezing point temperature is - 10°C and the mean tempera- 
lure -2 to -3°C 

Fi g u re 6 Futssan ce calculee ale vapor at eur . I m tempera t ure 
du point de congelation est de -KFC ei la temperature 
moyenne de —2a -S*C 



taken into the count in Equation (16). The fluid 
temperature increases in the purrip (heat flow to the 
loop is omitted): 



c,m. 



(17) 



The known input parameters are the evaporation 
temperature T c the mass flow rate of fluid m fy the air 
inlet temperature T,|, the mass flow rate of air m 4 . 
The surface areas A a and A c can be calculated from 
the reference state values. The refrigerant and air side 
heal transfer coefficients are assumed to be constant. 
The Equations (7) were solved using the standard 
Newton- Raphson iteration. 

Sample results for turbulent flow are shown in 
Figures 6 and 7. Figure 6 shows the capacity of the 
evaporator vs the mass flow rate. Figure 7 shows the 
changes of the pumping power. At a specified mass 
flow rate, the refrigerant does not have a very great 
influence on the capacity of the evaporator. Better 
heat transfer of PF is compensated by the lower 
specific heal, which gives a higher temperature change 




0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 
MASS FLOW RATE, kg/s 



1.2 1.3 



Figure 7 Calculated pumping power demand for the cases 
of Figure 6 

Figure 7 Demande de puissance de pompage calculee pour 
les cas present es dans la Figure 6 




Figure 8 Schematic of the experimental set-up. A, Storage 
tank: C. evaporator; D, condenser; E, plate heat exchanger; 
F, electric healer; G, plate finned air cooling coil 

Figure 8 Schema du dispositij experimental. A . Reservoir de 
stock age; C. ewporateur; D, condenseur; E. echangeur de 
chaleur a plaques; F. chauffage electrique; G. serpent in de 
re/roldissement d'uir a plauues-ailettes 
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ol* the refrigerant. However, if the mass How rate is 
chosen to gel the same pumping power, a higher 
capacity of the evaporator is obtained lor PI". For 
fluid mixtures having a lower freezing temperatures, 
the improvement is even more pronounced. 

Experimental comparisons 

Experiments were carried out in a test rig illustrated in 
Figure 8. The liquid -liquid heat exchanger was a 
standard plate type. The flow rale in both sides were 
adjusted at the same value in order to determine the 
heat transfer coefficient between plates and fluid. The 
air cooler was a plate fin heat exchanger having 5 mm 
fin spacing, 0.25mm fin thickness, I4.5/I3.5mm 
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Figure 9 Measured thermal conductance of a plate heat 
exchanger. The freezing poini temperature is - I0°C 
Figure 9 Conductance thermique mesuree dun echangeur de 
chateur a plaques. La temperature du point de congelation est 
de - IPC 



copper tubes with 36/30 mm staggered spacing. The 
air How rate was maintained constant. 

From the measured temperatures. How rates and 
pressures, pumping power demand and thermal 
conductance at various flow rates were determined. 
Some examples of the results are shown in Figures 9 
and 10, The temperature of the secondary refrigerant 
in the figures is -5 to 3°C. 

Figure 9 shows* the thermal conductance, the 
temperature change and the pumping power for the 
plate heat exchanger when the freezing temperature is 
- 10'C and the mean temperature -4 C. The thermal 
conductance G is defined by Equation (10). 

Figures 9 and 10 enable comparisons of the 
performance results in different ways. One way is to 
compare heal transfer and conductance values al the 
same pumping power as above in theoretical studies. 
For example, if we lake the value 20 for the pumping 
power, from Figure 10 the mass flow rate aboul 0.62 is 
obtained for EA and PG, 0.65 for EG and 0.72 for 
PF. Usinv jicse (low rales, ihe inductances are 
obtained using Figure 9. The conductance value is the 
lowest for PG. EA and EG give aboul 5% higher and 
PF about 60% higher conductances. In practice, the 
improvement is smaller because the secondary refrig- 
erant flows only in one side of Ihe heal exchanger. 
However, also in this case PF clearly improves heat 
transfer. 

Similar but smaller improvements were observed in 
the air cooling coil. This is due to low heat transfer at 
the air side and thus the weaker influence of the 
secondary refrigerant on the total resistance. In 
practice, often more important is the possibility to 
avoid the laminar flow in the tubes of the coil. 

Summary and conclusions 

The properties of the non-traditional secondary 
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refrigerant, PF, were studied theoretically and experi- 
mentally. The most important advantages compared 
to the traditional aqueous solutions of alcohols and 
glycolcs arc good thermodynamic properties, a low 
toxicity and non-ftammability. Volumetric heat is 
lower, which requires a higher rate of mass flow, if the 
temperature change of the fluid is maintained 
constant. However, because of clearly lower values 
of viscosity, the pumping power demand is at same 
level. Also, good thermal conductivity promotes good 
heat transfer. 

According to the results of theoretical and experi- 
mental comparisons, higher values of the thermal 
conductance of heat exchangers are achieved. One 
important feature in many cases is the possibility to 
avoid the laminar flow region in heal exchangers at 
lower velocities. 
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Appendix: the reference state 

r a2i0 =<rc 

7 co = -5"C 
m iQ =0.25kgs - ' 
rii tQ = 0.088 kgs"' 
<* a0 = 20Wkm" 2 

= 500 W km -2 
R u =0.05 km 2 W - ' 
R r = 0.002 km 2 W' 1 



